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^ I Spectral and photometric observations of nearby galaxies show a correla- 

I tion between the strength of their mid-IR aromatic features, attributed to PAH 

^ ' molecules, and their metal abundance, leading to a deficiency of these features in 

low-met alii city galaxies. In this paper, we suggest that the observed correlation 
^ ' represents a trend of PAH abundance with galactic age, reflecting the delayed 

I injection of carbon dust into the ISM by AGB stars in the final post-AGB phase 

of their evolution. AGB stars are the primary sources of PAHs and carbon dust 
in galaxies, and recycle their ejecta back to the interstellar medium only after a 
few hundred million years of evolution on the main sequence. In contrast, more 
massive stars that explode as Type II supernovae inject their metals and dust 
almost instantaneously after their formation. We first determined the PAH abun- 
dance in galaxies by constructing detailed models of UV-to-radio SED of galaxies 
that estimate the contribution of dust in PAH-free Hll regions, and PAHs and 
dust from photodissociation regions, to the IR emission. All model components: 
the galaxies' stellar content, properties of their Hll regions, and their ionizing 
and non-ionizing radiation fields and dust abundances, are constrained by their 
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observed multiwavelength spectrum. After determining the PAH and dust abun- 
dances in 35 nearby galaxies using our SED model, we use a chemical evolution 
model to show that the delayed injection of carbon dust by AGB stars provides 
a natural explanation to the dependence of the PAH content in galaxies with 
metallicity. We also show that larger dust particles giving rise to the far-IR 
emission follow a distinct evolutionary trend closely related to the injection of 
dust by massive stars into the ISM. 

Subject headings: ISM: dust - infrared: galaxies - galaxies: starburst - galaxies: 
evolution - stars: post-AGB - supernovae remnants 



INTRODUCTION 



Spectral and photometric observations of nearby galaxies with the Infrared Space Obser- 
vatory (ISO) and the Spitzer Space Telescope have provided the opportunity to investigate 
the inter-relations between global galactic properties, such as morphology, star formation 
activity, spectral energy distribution (SED), metallicity, and dust abundance and compo- 
sition. In particular, these observations have enabled detailed studies of the correlation of 
dust abundances and composition with metal enrichment in galaxies spanning a wide range 
of metallicities. Since the metallicity of galaxies evolves monotonically with time, galaxies 
with different metallicities provide snapshots of the evolutionary history of galaxies. 

An exciting result provided by ISO spectral observations of nearby galaxies was the 
discovery of a striking correlation between the strength of their mid-IR aromatic features 
and their metallicity flMadden et al.ll2006l ). Low- metallicity galaxies exhibited very weak 
or no aromatic features. Observations of the 8-to-24 /xm bands flux ratio obtained with 
Spitzer/IRAC and Spitzer/M IPS instruments showed a correlation of this flux ratio with the 
galaxies' oxygen abundance (lEngelbracht et al.ll2005l ). Since the Spitzer /IRACs^m band is 



supposed to trace the strength of the aromatic features, and the Spitzer/MIPS24^m that of 
the continuum emission from the hot non- aromatic dust component, this correlation seemed 
to confirin the t rends discovered by ISO. Thi s correlation was independently confirmed by 
Wu et al.l (120061 ) and lO'Halloran et al.l (120061 ) from spectral observations of low-metallicity 
blue compact dwarf galaxies with the Spitzer/ IKS instrument. 

The aromatic features are most c ommonly attributed to the vibrational rnodes of Poly- 
cyclic Aromatic Hydrocarbons (PAHs: lLeger fc Pugetlll984l : lAUamandola et al.lll985l ). which 
are large planar molecules made of 50 to 1000 carbon atoms. Their ubiquity makes them 



an important component of dust models (jPesert et al.lll990l : iDwek et al.l 119971 : IZubko et al. 
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2004: iDraine fc Lil 120071) . locking up ~ 15 — 20% of the total amount of interstellar carbon 
( iZubko et al.ll2004j . with solar abundance constraints). They mostly reside in galactic pho- 
todissociation regions (PDRs), wh ere they play an important r ole in the heating of th e gas 



by providing photo-electrons (e.g. iTielens fc HoUenbachl Il985l : iBakes fc Tielenj Il994l ). and 



in interstellar chemistry by providing surfaces for chemical reactions. Because of their small 
sizes, PAHs are stochastically heated by the interstellar radiation fields. The relative strength 
of some of the aromatic features depends on their ionizatio n state, and varie s therefore sig- 



nificantly with the physic al condition of the environment (IHony et al.l l200ll : IVermeij et al. 



2OO2I : iGalliano et al.l 120071 ). Understanding the evolution of PAHs and their relation to the 
global properties of galaxies is therefore extremely important because of the complex mutual 
influences between PAHs, and their ambient radiative and gaseous surroundings. 

Several explanations have been offered for the correlation of the intensity of the PAH fea- 
tures with metallicity. The first suggests that the tr end reflects an increase in t he destruction 



efficie ncy of PAHs in low metallicity environments (IGalliano et al.ll2003l . l2005l : iMadden et al. 



20061 ). Low-metallicity environments are bathed with harder photons than our Galaxy, due 
to the higher effective temperature of their stars, and their young age. Moreover, the paucity 
of dust allows this hard radiation field to penetrate deeper into the ISM, compared to high- 
metallicity systems, selectively destroying the PAH molecules by photoevaporation or pho- 
todissociation. T his explanat i on is consistent with the models of population synthesis and 
dust evolution of iDwek et al.l (120001 ). Assuming that PAHs are efficiently destroyed by UV 
photons in Hll regions, their models showed an evolutionary trend of PAH features with 
time, as the relative contribution of ionizing OB s tars to the galaxy ' s SEP decreases with 
time. A second explanation has been proposed by lO'Halloran et al.l (120061 ). who suggested 
that PAHs are destroyed by the numerous shocks observed in low metallicity systems. To 
support their proposition, they showed an anti-correlation between the PAH-to-continuum 
ratio and the [Fe ll]25.99;im/[Ne Il]i2.8i/im line ratio, the latter being supposedly a shock tracer. 
The problem with this explanation is that there is no observational ev idence that PAHs are 
selectively destroyed in shocks. On the contrary, iReach et al.l (120021 ) showed that, in the 
shocked medium of 30 391, both the PAH features and the underlying continuum disappear. 

All previous explanations attribute the paucity of PAHs to de structive pro cesses that 
are more efficient in the early stages of galaxy evolution. In contrast, iDwekl (120051 ) suggested 
that the observed correlation reflects an evolutionary trend of the sources of interstellar PAHs 
with metallicity. PAHs and carbon dust are mostly produced in asymptotic giant branch 
(AGB) stars which, unlike massive stars, recycle their ejecta into the ISM after a significantly 
longer time of main sequence evolution. The observed correlation of PAH line intensities with 
metallicity is therefore a trend of PAH abundance with galactic age, reflecting the delayed 
injection of PAHs and carbon dust into the ISM by AGB stars in their final, post-AGB, 
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phase of their evolution. 



Such distin c t evo lutionary trends of SN and AGB produced dust w ith time were pre- 
dicted by iDwekl (119981 ) and more recently by lMorgan fc Edmunds! (120031 ) . In particular, the 
latter show the evolutionary trend of AGB dust with time. This trend can be translated 
as a trend with metallicity if galaxies approximately share a global "cosmic" star formation 
history. 

The goal of this paper is to present a detailed evolutionary model to examine whether the 
observed trend of PAH line intensity with metallicity reflects an evolutionary trend of PAH 
abundance with metallicity. This requires the determination of PAH abundance from the 
strength of their aromatic features in the galaxies for which this trend has been observed, and 
the use of a chemical ev olution model to follow the change of PAH abundance with galactic 
metalhcity JDweklll998h . The paper is organized as follows. We first present in §2] the sample 
of nearby galaxies that were considered in our analysis. In ^ we describe the method we 
used to separate the contributions of H ll regions and photodissociation regions to the global 
SED, in order to determine the abundances of PAHs and larger grains in these galaxies in 
^ In ^ we briefly describe the chemical evolution model used in the calculations, and 
compare its results to the abundances derived from our SED modeling. The results of the 
paper are briefly summarized in ^ 



Throughout this paper we will refer to the solar abundances by iGrevesse fc Sauval 

( 119981 ). the oxygen number abundance being 12 + log(0/H)Q = 8.83, the Helium and heavy 
elements to gas mass ratios Yq = 0.248 and Zq = 0.017, respectively. Besides, we assume 
that the helium abundance is independent of the metallicity. 



THE SAMPLE OF NEARBY GALAXIES 



2.1. Source Selection 



In order to properly characterize the PAH emission, we considered galaxies whose 
mid-IR spectrum has been observed, either with one of the spectrographs onboard the 
ISO satellite, or with the Spit zer/IRS. We combin ed ISO sam ples of starbursts and AGNs 
( iLaurent et al.l 120001) . spirals (IRoussel et al.l l200ll ) , ellipticals (IXilouris et al.l |2004| ) , dwarf 
galaxies (IMadden et al.1 120061), and the low-metallicity sources of the Spitzer sample pre- 
sented by lEngelbracht et al.l (120051 ) which were observed by Spitzer/IKS. Among the ellip- 



ticals, only NGC 1399 has a mid-IR spectrum. Consequently, our sample includes various 
types of galaxies and covers a wide range of metallicities and star formation activity (Ta- 
ble [I]). 
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The modeling that will be presented in ^ requires the assembly of data covering the 
stellar as well as the dust emission components, for each galaxy. First, we need most of 
the U, B, V, R, I, J, H, K fluxes, to constrain the stellar spectrum shape. Photometric 
J, H and K bands are available for almost all of our sources, thanks to the 2MASS survey 



fiJarrett et al.ll2003D . We rejected the galaxies Arp 118, Arp 236, Arp 299 and NGC 4038, 
for which no B and V bands were reported. Second, the far-IR SED is used to constrain the 
intensity of the interstellar radiation field (ISRF). Hence, we rejected galaxies which have 
not been detected by IRAS or MIPS, such as HS 0822+3542, Tol 1214-277 and Tol 65. In 
addition, our mass estimates are normalised to the Hydrogen mass. We therefore rejected 
the galaxy IRAS 23128-5919, for which no H I observation was reported. Finally, we removed 
M 31 from our sample, since its angular diameter is too large to build a consistent observed 
total SED. 

The global properties of the selected sources are presented in Table [1] If relevant, the 
distances were homogenised to Hq = 71 kms~^ Mpc~^. The masses have been scaled to the 
adopted distance. 



Table 1. Select Properties of the Sample 



Name 


R.A. 


Dec. 


Mid-IR 


Distance 


12 -f 


log(0/H) 


Mh, 




Mr, 




Notes 




(J2GGG) 


(J2GGG) 


spectrograph 


(Mpc) 




[ref.] 


(1G8 Mq) 


[ref.] 


(108 Mq) 


[ref.] 




Haro 11 


00''36'"52!5 


-33°33'19" 


Spitzer/IRS 


92 


7.9 


[2] 


< iCO 





^ 1 





Pec H II 


NGC 253 


00'H7™32!9 


-25°17'18" 


ISO/CAM 


3.3 


9.0 


[3] 


18 


[4] 


17 


[5] 


Hii 


NGC 520 


0l''24™34'?9 


-fG3°47'31" 


ISO/CAM 


27 


(•) 




35 


[6] 


35 


[7] 


Pec H II 


NGC 613 


0l''34™17;5 


-29°24'58" 


ISO/CAM 


19 


9.2 


[8] 


37 


[9] 






Sy 


NGC 891 


02''22™33!4 


-h42°20'57" 


ISO/CAM 


9.6 


8.9 


[10] 


76 


M 


49 





Edge-on 


NGC 1068 


02''42™40';6 


-GG°00'47" 


ISO/CAM 


15 


9.0 


[11] 


22 


[12] 


74 


[13] 


Sy 


NGC 1097 


02''46™19n 


-3G°16'28" 


ISO/CAM 


12 


9.0 


[14] 


51 


[15] 


^ 7.3 


[16] 


Sy 


NGC 1140 


02''54™33'?5 


-1G°01'44" 


ISO/CAM 


25 


8.0 


[17] 


52 


[18] 






Irr H II 


NGC 1365 


03''33™35';6 


-36°08'23" 


ISO/CAM 


19 


9.1 


[19] 


13G 


[20] 


170 


[21] 


Sy 


SBS 0335-052 


03''37™44'?G 


-G5°02'38" 


Spitzer/IRS 


54 


7.3 


[22] 


9.9 


[23] 






BCD 


NGC 1399 


03''38™29n 


-35°27'03" 


ISO/CAM 


21 


CO 




< 2 


M 






cD 


IC 342 


03''46™49!7 


-l-68°05'45" 


ISO/CAM 


3.8 


8.9 


[24] 


130 


M 


50 





Hii 


NGC 1569 


04''30™49!1 


-F64°50'52" 


ISO/CAM 


2.2 


8.2 


[25] 


1.3 


[26] 


0.50 


[27] 


Irr H II 


NGC 1808 


05''07™42';3 


-37°30'47" 


ISO/CAM 


11 


9.1 


[28] 


18 


[29] 


20 


[30] 


Sy 


II Zw 40 


05''55'"42';7 


-fG3°23'29" 


ISO/CAM 


10 


8.1 


[31] 


4.2 


M 


0.23 





BCD 


He 2-10 


08''36™15'?2 


-26°24'34" 


Spitzer/IRS 


8.7 


8.9 




3.1 


[32] 


1.6 


m 


Irr H II 


I Zw 18 


09''34™02';G 


-|-55°14'28" 


Spitzer/IRS 


13 


7.2 




1.2 


[33] 






BCD 


M 82 


09''55™51'?8 


-f69°40'46" 


ISO/CAM 


3.6 


9.0 


[34] 


9.0 


[35] 


13 


[36] 


Irr H II 


NGC 3256 


10''27™5in 


-43°54'17" 


ISO/CAM 


37 


8.9 


[37] 


62 


[38] 


300 


[39] 


Pec H II 


Mrk 33 


10''32'"31':9 


-h54°24'04" 


Spitzer/IRS 


20 


8.4 


m 


4.3 


[40] 


0.70 


[41] 


Irr H II 


Mrk 153 


10''49™05';G 


-H52°20'08" 


Spitzer/IRS 


37 


7.8 


[42] 


6.7 


[43] 






BCD 


VII Zw 403 


ll''27"59!9 


+78°59'39" 


Spitzer/IRS 


4.5 


7.7 


[22] 


0.69 


m 






BCD 


UM 448 


ll'H2™12M 


-fGG°20'03" 


Spitzer/IRS 


70 


8.0 


[44] 


47 


[45] 


24 


m 


Pec H II 


NGC 4945 


13''05™26';2 


-49°28'15" 


ISO/CAM 


3.9 


(?) 




45 


[46] 


17 


[47] 


Edge-on Sy 


Centaurus A 


13''25™28^G 


-43°01'06" 


ISO/CAM 


3.8 


- 9 


[48] 


11 


[49] 


2.0 


[50] 


Elliptical Sy 


M51 


13''29™52';7 


-H47°11'43" 


ISO/CAM 


8.4 


8.7 


[51] 


50 


[52] 


56 


[53] 




M83 


13''37™00!7 


-29°51'58" 


ISO/CAM 


4.5 


9.2 


[54] 


51 


[55] 


49 


[56] 


Hii 



Table 1 — Continued 



Name 


R.A. 


Dec. 


Mid-IR 


Distance 


12 -F 


- log(0/H) 


Mhi 




Mr, 




Notes 




(J2000) 


(J2000) 


spectrograph 


(Mpc) 




[ref.] 


(10» Mo) 


[ref.] 


(10« Mq ) 


[ref.] 




Tol 89 


14''01™21';5 


-33°03'50" 




15 


8.0 


[57] 


13 


[58] 






Hn 


Circinus 


14/113^09^6 


--65°20'21" 


ISO/CAM 


4.0 


(?) 




55 


m 


11 


[59] 


Sy 


NGC 5253 


13''39™55!7 


-31°38'29" 


ISO/SWS 


3.3 


8.2 


m 


0.91 


[60] 


<o.i 


[61] 


Irr H II 


Arp 220 


15''34'"57;2 


+23°30'11" 


ISO/CAM 


73 


(?) 




280 


m 


160 


[62] 


ULIRG 


NGC 6240 


16''52™58':8 


+02°24'06" 


ISO/CAM 


98 


(?) 




88 


m 


310 


m 


LIRG 


NGC 6946 


20''34™51^2 


+60°09'17" 


ISO/CAM 


5.5 


9.1 


[63] 


56 


[64] 


33 


[65] 


Hii 


Mrk 930 


23''31™58';3 


+28°56'50" 


Spitzer/IRS 


73 


8.1 


m 


28 


[66] 






Hii 


NGC 7714 


23''36™14!1 


+02°09'19" 


Spitzer/IRS 


37 


8.5 


[67] 


17 


[68] 


22 


m 


Pec H II 
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2.2. ISO/CAM Data Reduction 



Most of the sources in Table [D were obse rved witli ISO /CAM (ICesarskv et al.l 119961) 



on bo ard tlie ISO satellite flKessler et al.l Il996l ) . These spectra were used by iMadden et al. 



(120061 ) . and we refer to this work for a detailed description of the data reduction. The CVF 



performed spectral imaging using a 32 x 32 detector array, with a sampling of 3"pixel~^ or 
6"pixel~^ in our cases, from A = 5 /im to 16.5 /im with one pointing of two CVFs, from 
A = 5 to 9.5 /im and from A = 9.0 to 16.5 /im. The spectral resolution goes from A/AA = 35 
to 51 across the full spectra. 

For the d ata treatment, we used the CAM Interactive Reduction (CIR, v ersion AUGOl; 



Chanialll2003l ). The subtraction of the dark currents was performed using the lBiviano et al. 



(119981 ) model which predicts the time evolution for each row of the detector, taking into 
account drifts along each o rbit and each revo lution. We masked the glitches using multi- 
resolution median filtering (jStarck et al.lll999l ) on each block of data after slicing the cube. 
Additional deglitching was performed manually, examining the temporal cut for each pixel. 
We c orrected the systematic memory effects using the Fouks- Schubert method (ICoulais &: Abergel 
2000l ). We computed a hybrid fiat-field image placing a mask on the source and computing 
a flat field outside this mask from the median of the temporal cut for each pixel. For the 
pixels which were on-source, the flat-field response was set to the corresponding calibration 
flat-field. The conversion from Analog Digital Units to mJy/pixel was performed using the 
standard in-flight calibration data base. To remove the sky contribution, the sources smaller 
than the array were masked and, for a given wavelength, the median of the pixels which are 
off-source were subtracted from each pixel. For the more extended sources, we subtracted 
an independently observed zodiacal spectrum. The intensity of this spectrum was a free 
parameter varied in order to match the properly sky subtracted fluxes in the LW2 (6.7 /im) 
and LW3 (14.3 /xm) broadbands. The final product is a 3D spectral-image of each galaxy. 
We integrated the spectrum into an aperture encompassing the entire galaxy, and to obtain 
the global SED of the galaxy. When the angular size of the source was larger than the one 
of the array, we scaled-up the spectrum to match the IRASu^m broadband flux. 



2.3. Spitzer/mS Spectrum Extraction 



Several of the low-metallicity sources in Table [T] were not observed by ISO/CAM, 
we therefore complemented our database with publically rele ased mid-IR spectra from the 



Spitz er/ IRS spectrometer on board the Spitzer Space Telescope (IHouck et al.ll2004bl : IWerner et al. 



2004h . Among these galaxies, the spectra o f SBS 0335-052, NGC 7714, VII Zw 403, Ha ro 11 
and I Zw 18 are described in details in iHouck et al.l (l2004al ). iBrandl et al.l (120041 ) and 
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Wu et al.l (120061 ) ■ We considered only low- resolution data, taken with the SL (Short-Low) 



module, from A = 5.2 fim to 14.5 /im, and the LL (Long-Low) module, from A = 14.0 fim to 
38.0 fim, both with a spectral resolution of A/AA ~ 64 — 128. 

We retrieved the Basic Calibrated Data (BCD) that have been preprocessed by the 
Spitzer Science Center (SSC) data reduction pipeline, and converted to flux density, corrected 
for stray light and flatfielded. The extraction of the spectra from the 2D space/wavelength 
images was p erformed with the S pectral Modeling, Analysis and Reduction Tool (SMART, 
version 5.5.6; iHigdon et al.ll2004l ). We first inspected the BCD images and identified the 
hot pixels which have not been masked out by the SSC. We replaced them by the median 
of their neighbors. For each module, an off-source position is observed. We subtracted this 
spectrum from the on-source one, in order to remove the sky emission. The extraction of 
the ID spectrum was performed inside a column whose width vary with the wavelength. 
We have excluded the bonus order. Then, the various frames, for each nod position, were 
coadded. Since the long wavelength end of the SL module and the short wavelength end of 
the LL one were not systematically overlapping, we finally scaled the SL module, in order to 
obtain a continuous spectrum. This scaling factor can be as large as 50%. At the time when 
this publication is written, the Spitzer/IKS data handbook recommends not to derive the 
signal-to-noise ratio from the uncertainties generated by the pipeline. Instead of that, we 
adopted the recommended systematic error of 20%. However, the SL module of the noisiest 
spectra (VII Zw 403, Mrk 153, Mrk 930, I Zw 18, Tol 89) exhibits fluctuations larger than 
this value. To take into account these statistical variations, we smoothed these spectra into 
a AA ~ 0.2 fim window (4 points), and took the standard deviation inside this window as 
the error. Similarly to what we did with ISO/ CAM spectra ( §2.2p . in order to compensate 
the fact that we may be overlooking some extended emission, we scaled-up the Spitzer/IKS 
SL and LL spectra to match the IRASufim and IRAS25^m- In the case o f II Zw 40, the 
Spitzer/IRS spectrum does not exhibit the PAH features that lMadden et al.l (120061 ) detected 
in the extended emission, thus we prefer to use the ISO /CAM data for this galaxy. For 
Tol 89, the Spitzer/lKACsum broadband flux is higher than the integrated spectrum into 
the same band. This is probably due to the fact that the slit of Spitzer/IRS measures only 
the nucleus emission, which likely has a steeper continuum and weaker features than the 
extended emission contributing to the total broadband. Therefore, we will consider the mass 
of PAHs derived for this galaxy to be a lower limit. 
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2.4. The Multiwavelength SEDs 



The UV-to-radio SED of each galaxy was buih using data from different catalogs. We 
extensively used the HYPERLEDA and NE D databases. Most o f the SEDs includ e IRAS 
broadband fluxes at 12, 25, 60 and 100 u m jMoshh fc et al.lll99ol : iRice et aLlEoSsf ). and J, 
H, K 2MASS fluxes f lJarrett et al.ll2003l ). For som e of the galaxies, we completed the IR 
SED with the broadband Spitzer data reported by iDale et al.l (120051 ). We kept only data 
encompas sing the entire galax y. The optical data were corrected for Galactic extinction, 
using the ISchlegel et al.l (119981 ) extinction maps. 



For NGC 5253, we used the ISO/SWS observations of ICrowther et al.l (119991 1 to model 
the spectrum of this galaxy. 



3. MODELING THE PANCHROMATIC SEDS OF GALAXIES TO 
DERIVE THE PAH AND DUST ABUNDANCES 



3.1. Motivations and Approach 

Determining the PAH abundance from the 8 and 24 fun fluxes or even from the mid-IR 
spectrum requires: (i) knowledge of the contribution of other dust species to the mid-IR 
emission, in order to determine the intrinsic emission in the PAH bands; and (ii) knowledge 
of the ISRF they are subjected to, since PAHs are stochastically heated. The dust that 
contributes to the mid-IR emission is hot dust radiating at equilibrium temperature in H ll 
regions, as well as very small grains undergoing temperature fluctuations in PDRs. Therefore, 
determining the PAH abundance requires a panchromatic approach, to model the SEDs of 
galaxies. It requires knowledge of star formation history, and processing of the stellar SED 
by the gas and dust in the ionized and neutral phases of the galaxy. 

The 8/24 fim flux ratio provides information on the strength of the PAH features in the 
galaxy, although it is relatively limited. In general, the Spitzer/IRACs^m. flux is dominated 
by the 6 — 9 fim PAH features, but in some cases, such as low metallicity galaxies or 
deeply embedded sources, the 8 p^m fl ux may be dominated by, respectiv ely, the continuurn 



(iMadden et al.l 



20061 : IWu et al.l 120061 ) or the 9.7 /im silicate feature (e.g. ISpoon et al.l 12004 



Hao et al.l 120071 ). When there are no PAHs and the Spitzer/IRACs^m band is dominated 
by the continuum, the 8/24 /im ratio reflects the color of the mid-IR emitting silicates and 
graphites. The Spitzer/M.IPS2Afj.m flux is dominated by the continuum emission from the 
silicate and carbon g rains. However , this part o f the SED is very sensitive to the abundance 
of small grains (e.g. iGalliano et al.l l2003l . 120051 1. and can originate additionally from large 
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grains close to strong sources of radiation (e.g. iPlante fc Sauvagd l2002l : IVanzi fc Sauvage 



2004h . 



Fig. [T] shows the ratio of the two bands, for the sources presented in Table [H The figure 
shows that when the PAH-to-continuum ratio is small, the 8/24 fim ratio plateaus around a 
value of ~ 0.1, which roughly represents the color of the mid-IR continuum. The figure illus- 
trates that the separation between the PAH bands and the continuum is important to study 
the evolution of their strengths with the metallicity. The error bars on the ratio, in Fig. [H 
come from the propagation of the observational errors. The error bars on the metallicity 
were not systematically given by the authors who published them. An error of 0.1 dex in 
the 0/H number abundance accounts for the typical dispersion between independent mea- 
surements. There were 6 sources, in Table [H for which no metallicity measurements were 
reported. However, all of them have the morphology, optical colors, IR emission and H2 con- 
tent of solar or supersolar systems, with NGC 1399 probably having the largest metallicity. 
We assign an arbitrary Z = (2±1) Zq to these galaxies, and consider them to be uncertain. 



Knowledge of the shape of the ISRF to determine the PAH abundance is important, 
since the absorption properties of PAHs are significantly different from the ones of solid-state 
grains (carbonaceous or silicates) believed to be the carriers of the mid-IR continuum. Fig. [2] 
shows the different wavelength-dependence of the mass absorption coefficient, ^ahs, of various 
types of grains. It shows that the value of /tabs for PAHs drops by 4 orders of magnitude 
between 0.1 and 1 fim, where the stars emit most of their energy, while the Kabs of graphite 
and silicate dust drop by less than 2 orders of magnitude. However, most of the energy is 
absorbed at shorter wavelengths. Thus, the PAHs are more sensitive to very young stellar 
populations than the grains responsible for the 24 fim continuum (Fig. [3]), and a correct 
determination of their excitation rate should take into account this property. 



3.2. The Contribution of Hll Regions to the SED of a Galaxy 



Models of the spatial distribution of the dust in our Galaxy i ndicate that Hll re- 
gions dominate the infrared SED up to wavelengths of ~ 60 fim (jSodroski et al.l 119971 : 



Paladini et al.l 120071 ). In principle, the contribution of Hll regions could be derived from 



any tracer of the recombination rate, like the radio free-free continuum. However, the free- 
free-to-IR r atio depends significantly on the density a nd dust-to-gas mass ratio in the H ll 



Petrosian et al. 



region (e.g. 
regions as a function of these parameters. 



19721 : iMezger et al.lll974l ). requiring detailed modeling of Hll 
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Fig. 1. — Mid-IR color as a function of the oxygen number abundance, for the sample in 
Table [H This figure illustrates the importance of knowing the mid-IR continuum emission 
in order to derive the intensity of the PAH bands. Our observed SEDs have been integrated 
over the Spitzer/lRACsfim and Spitzer/MlPS24:^m bandpasses, in order to produce this plot. 
The two open circles are the (U)LIRGs of our sample. The grey stripe is the ztlcr linear 
correlation between log(0/H) and log(-F,^(8 fim) / F^{24 /im)). The hatched stripe shows the 
range of the ratio (0.06 < F^(8 yum)/Fj,(24 /im) < 0.13), which is a measure of the color of 
the silicate and graphite continuum, when the PAH features are weak. 
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Fig. 2. — The wa velength-de pendence of the mas s absorption c oefficient of the PAHs 
(IDraine fc Lil 120071 ) , graph ites flLaor fc Draind Il993l ) and sihcates f lWeingartner fc Draine 



200ll ). integrated over the IZubko et aLl (120041 ) size distribution, for the bare grain model 



with solar abundance constraints (BARE-GR-S). 
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Fig. 3. — Sensitivity of the mid-IR SED to the age of the stellar populatio ns. The two 



SEDs come from the same grain populations (abundances and size distribution; IZubko et al. 



20041 ): the only difference is the shape and intensity of the interstellar radiation field they 
are experiencing. The black SED is heated by a young UV-rich instantaneous burst, while 
the grey one is the subjected to the local ISRF. It shows that for a given far-IR spectrum, 
the PAH emission can vary by a factor of ~ 3 depending on the shape of the ISRF, therefore 
affecting the estimate of the PAH abundance by the same amount. 
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We consider H ll regions to be made of a central ionizing star, surrounded by PAH-free 
dusty ionized gas. The output of an individual Hll region therefore consists of: (i) escaping 
non- ionizing stellar radiation; (ii) IR emission from the reprocessing of the stellar emission 
by dust; and (Hi) radio free-free emission from the reprocessing of ionizing photons by the 
gas. We created a library of H ll region spectra for stars of different spectral types, embedded 
into a gas of different densities and dust-to-gas mass ratios. The total contribution of H ll 
regions to a galaxy's SED is derived by integrating the spectra of indvidual H ll regions over 
the stellar initial mass function (IMF). 

There are several simplifying assumptions in our procedure. Individual H ll regions are 
considered to be spherical, homogeneous, consisting of pure Hydrogen. In calculating the 
infrared emission, we neglected scattering of stellar radiation by the dust. We assumed 
that all the dust radiates at equilibrium dust temperature, and that PAH are destroyed and 
therefore absent into the Hll region itself. In calculating the emission from the ensemble 
of Hll regions, we assumed that they are not clustered, so that the total radiation is an 
IMF weighted sum over all spectral types, assumed to be surrounded by gas with the same 
density. 

The theoretical evolution of the SED of Hll regions and t h eir surrounding mole cular 
material has been extensiv ely discussed by iBressan et al.l fl2002l ) , iPanuzzo et al.l (120031 ) and 
Dopita et al.l (l2005l . l2006a| ]bl). In these studies, the spectrum of an Hll region depends on 
its evolutionary state, controlled by its age, density and the dissipation time scale of the 
surrounding molecular clouds. In comparison, our approach is more flexible and empirical, 
in which the parameters can vary more freely instead of being constrained by evolutionary 
models. 



In the following, we describe how we derive the SED of H ll regions. The definition of 
each mathematical variable is summarised in Table [21 



Table 2. Mathematical Variables Used Throughout the Paper. 



Name 


Units 


Definition 





General Variables 



v Hz Frequency of the radiation 

A /xm Wavelength of the radiation 

Mgas Mq Gas mass (Hydrogen, Helium and metals) of the region or galaxy 

M(just Mq Dust mass of the region or galaxy 

Zdust = -^^dustA^gas Dust-to-gas mass ratio 

Zgas Metal mass fraction in the gas phase (or metallicity) 

TTiH = 1-67 X 10"^"^ g Mass of a single H atom 

TTidust g Average mass of a dust grain 

a jjim Grain radius 

/(a) grain /im~^ Grain size distribution, normalized to unity 

Qabs(fi) ^) Absorption efficiency of a grain of radius a 

/•oo 

Cdust(i^) = / 7r'2^(5abs(a, z^) /(a) da cm^grain"^ Average dust cross-section per grain 
Jo 

m Mq Mass of individual stars 

(Dim) Mq^ Initial mass function 



H II Region Variables ( ^^21) 



ui = 3.29 X 10 Hz Frequency of the Lyman break 

Rs cm Radius of the equivalent dust-free Stromgren sphere of a single H ii region 

y = r/Rg Radial coordinate normalised to Rs 

TT-H H cm~^ Number of H atoms (or protons) per unit volume 

ndust grain cm"'^ Number of dust particles per unit volume 



Table 2 — Continued 



Name 


Units 


Definition 






e cm 


Number of electrons per unit volume 




x{y) = He/nn 




Ionization fraction of Hydrogen 






2 

cm H 


Ionization cross-section of Hvdroeen fEo. 2.4 of Soitzer 1978) 




T' 




Electron temperature 




rp 


is. 


Equilibrium dust temperature 






Lq Hz 


iotal monocliromatic luminosity, passing through the shell of radius y 






Lq Hz-1 


Intrinsic monochromatic luminosity emitted bv the central star, given bv Panagia ( 


1973) 


LBiv) 




Luminosity of the case B recombination lines at a radius y 






Lq Hz-1 


Monchromatic luminosity reprocessed by the dust at a radius y 








Monchromatic luminositv of the free- free cooling at a radius v fEa. 3.54 of Soitzer 


1978) 


aA{Te) 


■? —1 

cm s 


Case A recombination coefficient fTable 2.1 of Osterbrock 1989) 






s —1 
cm s 


Case B recombination coefficient fTable 2.1 of Osterbrock 1989) 




L^^\i^,m,ZdvLst,nu) 


Lq Hz-1 


Monochromatic luminosity of an individual H ii region 




^burt 


Myr 


Age of the burst of star formation 




L^^^^^lu, tbursti -^dusti ^h) 


Lq Hz-1 


Monocromatic luminosity of a distribution of H II regions 




Global SED Variables (g33D 




Myr 


Age of the galaxy 






Mq 


Mass of the non-ionizing stars 




/sync 




Ratio between the synchrotron and the free-free at A = 1 cm 






Mq 


Dust mass in the PDRs 




/pah = AfpAH/Mdust 




PAH mass fraction 




/+ = MpAH+/MpAH 




Fraction of ionized PAHs 




Li 


Le 


Luminosity of a mid-IR ionic line 





Table 2 — Continued 



Name 



Units 



Definition 



U = / Uy{v) diy 
Jo 

L*iu, tsF,M^) 

Lj/ (z/, /sync) 
L^-b(j,,[/,Md,st,/pAH,/+) 

LPDR(z., a, [/„,[/+, Mdust, /pah, /4 

-^'burst ~ 



(-Pe-^burst) 
(-fe-Z^PDR) 



pH unburst 


lPDI^ di. 

pPDR^PDR 




L* dz^ 



ergs ""^ 


cm Hz 


erg s ^ 


cm 


Lq Hz- 


-1 




-1 


i-/0 xlZ 


-1 


Lq Hz- 


-1 


Lo 








Lo 




Lq 




Lq 




Lq 





Monochromatic radiation density in the PDRs 

Integrated radiation density in the PDRs 

Monochromatic luminosity emitted by the non-ionizing stars 
Monochromatic luminosity of the synchrotron radiation 
Monochromatic luminosity emitted by the dust in a given U 
Monochromatic luminosity emitted by the dust in PDRs 
Escaping fraction of the H ii region photons 
Escaping fraction of the PDR photons 

Intrinsic luminosity radiated by the H ii regions 
Escaping luminosity from the Hii regions 
Luminosity radiated by the dust in the PDRs 
Escaping luminosity from the PDRs 
Intrinsic luminosity radiated by the non-ionizing stars 
Escaping luminosity from the non-ionizing stars 



Elemental and Dust Evolution Variables 



^gas it) 
^gas 
5^SFR(i) 
■^ISM 

r(m) 

mej (t) 

Yz{m) 

/^gas(*) = Sgas(i)/Sg3 



M0PC-2 

Mqpc-2 
Mq yr~^ pc~^ 



Myr 

Mq 
Mq 



Gas mass surface density 

Initial gas mass surface density 

Star formation rate surface density 

Total metal mass fraction (gas and dust) in the ISM 

Lifetime of a star of mass m 

Mass of gas ejected by a star of mass m, after a time r(m) 
Mass of metals ejected by a star of mass m, after a time r(m) 
Reduced gas mass of the system 
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3.2.1. The Modeling of Individual Hu Regions 

We assume that PAHs are totally depleted inside the ionized gas phase (the actual 
Hll region), and that their mass fraction is constant outside of these regions, inside the 
photodissociation regions. Furthermore, we assume that the Hll regions are homogeneous, 
spherica l, and that they contain only H ydrogen and dust. The dust-free Stromgren radius 
Rs (e.g. ISpitzerlll978l : lOsterbrockl Il989l ) of such a region is defined by the balance between 
the rate of ionizing photons emitted by the central star and by the electrons recombining to 
the ground level, and the rate of recombinations to any level higher than the fundamental 
state: 

^Rlnla^{T,) = l^^^^du. (1) 

We use this quantity to normalize the radial coordinate r to the dimensionless radius y = 
r/Rg. In each shell of radius y and thickness dy, the Hydrogen optical depth is given by: 

dTu{i^, y) = [l- x{y)] Rs x an{T^) dy, (2) 

the dust optical depth is: 

drdust(i^, y) = '^dust(y) Rs X o-dust('^) dy, (3) 

and the volume of the shell is dV{y) = An y'^ dy. We solve the radiative transfer equation, 
from the central star to the photoionization front: 

dL^{u, y) _ f drH(^^,y) ^ d7dustOAy)\ ^ dL^"'^^(z/, ^ dL^Jji^, y) 



dy \ dy dy J dy dy 

together with the photoionization equilibrium in each shell: 

aA{Te) x{yfn^ dV{y) x hui = / L^{u, y) drH(i^, y) du. (5) 

J Vl 

Eq. ([5]) gives the value of x{y). 

We assume that all the case B recombination lines are resonantly scattered by the gas 
and are finally absorbed locally by the dust. Therefore the dust is heated both by the stellar 
continuum and these recombination lines: 



/•oo 
JO 



drdust(i^, y) LB{y) = ndust^(y) / crdust(«^) 47r5j,(i/, Tdust) du. (6) 

dy Jo 



dy 



In Eq. ([6]), we assume that the dust is at thermal equilibrium with the radiation field. Very 
small grains may still undergo temperature fluctuations, but since their size distribution is 



Table 2 — Continued 



Name 


Units 


Definition 


Tdustit) 

{miSM) 


Myr 

Mo 


Dust lifetime 

Average gas mass swept-up by a single SN il 
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poorly known, we neglect this effect in calculating the IR emission. This will lead to an 
underestimate of the ~ 1 — 20 fim continuum. We will discuss in §3.31 an empirical way to 
compensate this underestimation. 

Finally, we t ake into account the fact that dust sublimates in the center of H ll regions 



[e.g. Ilnouell2002l ). by considering that each shell where the equilibrium temperature of the 
dust exceeds its sublimation temperature is free of this dust specie. We a dopt sublimat ion 
temperatures of 2500 K and 1800 K for graphite and silicate, respectively (jKruegelll2003l ). 



We solved Eqs. and ([5]) for stars of masses 8 Mq < m < 100 Mq, with a grid of 
densities ranging from tt-h = 1 cm~'^ to rin = 3 x 10^ cm~^, and a g rid of dust-to-gas mass 



ratios ranging from Z^ust = 1/120 (Galactic value; IZubko et al.ll2004l ) to ^dust = 1/12000 



3.2.2. The Total Contribution of Hu Regions 

At the scale of a galaxy, the SED of the ionized gas phase is the combination of several 
H II regions. We assume that all the H ll regions of a given galaxy have the same density, 
and the same dust-to-gas mass ratio. We adopt a Salpeter initial mass function: 

0(m) oc m~^'"^^ for m_ < m < m+ (7) 

normalised to / 0(m) dm = 1 



with 



m_ = 0.1 Mq 
m+ = 100 Mq 



The total SED of an ensemble of H II regions is then given by: 

/•'^+(*burst) 

L^"'*(i^, tburst, Zdnst, nu) = iTi, Zdust, ^h) X 0(m) dm. 



The upper mass decreases with the age of the burst of star formation. Fig. H] shows several 
of these SEDs. 

In our case, the age of the H ll region does not have an important effect on the shape 
of the Hll SED. Indeed, the age affects essentially the ionizing-to-non-ionizing photon rate. 
Since our H ll SEDs are integrated only out to the photoionization front, the ionizing photons 
dominate the total power input. Therefore, in what follows, we will consider Hll SEDs with 
an age of 4 Myr. 
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1 100 10* 10® 

Wavelength [jum] 

Fig. 4. — SED of an ensemble of Hll regions (Eq. [8]): for a given age and dust-to-gas mass 
ratio, varying the density (top panel); and for a given age and density, varying the dust-to- 
gas mass ratio (bottom panel). The density tends to increase the dust temperature, as well 
as the free-free-to-IR ratio (top panel). Indeed, when the density is high, the dust absorbs 
directly a significant fraction of the ionizing photons. Inversely, when Z^nst drops, more 
ionizing photons are absorbed by the gas. The free-free-to-IR ratio is then higher, and the 
dust luminosity lower (bottom panel). 



3.3. The Inclusion of Emission from PDRs and the Resulting Galactic SEDs 

3.3.1. The Building Blocks 




1 100 10* 10' 

Wavelength [ju.m] 



Fig. 5. — Demonstration of the panchromatic SED model applied on a galaxy. 

A significant fraction of the radiation escaping from H ll regions is absorbed by the dust 
in the surrounding molecular clouds, as well as in the diffuse ISM. These dense and diffuse 
PDR components are also heated by the ambient non-ionizing radiation field. Technically, 
all the parameters of our model are adjusted simultaneously. However, for clarity, we can 
decompose the procedure in the following steps: 



1. Observations of at least two data points sampling the radio continuum constrain the 
synchrotron-to- free-free ratio. The synchrotron emission is L^^"'^ oc z/~°'^. 

2. The resulting free-free continuum, together with observations of the mid-lR continuum 
between ~ 5 and ~ 60 /im are used to constrain the parameters of the H ll regions. We 
use the Hll templates presented in §3.21 attenuated with a slab extinction P^^^{v) = 
exp [— rdust(i^)], to account for absorption by intervening dust to the Hll region. The 
power absorbed, Lburst — (-Pe-^^burst), contributes to the heating of the dust in PDRs. 
As mentioned in §3.2[ we have not taken into account the emission from stochastically 
heated grains in this phase. To correct empirically for this effect, we redistribute a 
part of the H II dust power into the sum of two modified black bodies of temperatures 
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130 K and 300 K (peaking around 15 and 5 /um). The luminosity of each individual 
black body is free to vary. We enforce that this redistribution does not exceed 40% of 
the total H ll dust luminosity. The dust-to-gas mass ratio in the H ll phase is assumed 
to be equal to the one in the PDRs. 

Optical/near-IR broadbands are used to constrain the esc aping diffuse stellar radia- 



tion, L*, using the stellar population synthesis code PEGASE (IFioc fc Rocca-Volmerange 



19971 ). with a Salpeter IMF, and a star formation rate as a functio n of time propor tional 



to a power-law of the gas surface density with an index of 1.4 (IKennicuttl Il998l ) . We 
vary the initial gas mass surface density T,^^ from 3 Mq pc~^ to 300 MqPc"^. Here 
again, we assume a slab extinction, for the sake of simplicity. The power absorbed, 
— (PgL^), contributes to the heating of the dust in PDRs. 

4. The far-IR/submm observed SED constrains the dust emission from the PDRs. In 
order to account for variations of the radiation density, U, in these regions, we assume 



a pow er-law distribution, following the empirical prescription introduced by lDale et al. 



fl200l[ ) 



L™(z/) oc / Ll''^{iy,U) X f/-° dU, (9) 
'u- 



where L^"'^(i^, U) is the dust SED corresponding to a single radi ation density U. W e 



adopt the dust properties of the Galactic diffuse ISM modeled by lZubko et al.l (120041 ). 
for bare grains with solar abundance constraints. We fix the silicate-to-graphite mass 
ratio, but let free to vary the PAH-to-dust mass ratio, /pah? as well as the fraction 
of ionized PAHs, /+. These PAH properties are constrained by the detailed fit of the 
features seen on the mid-IR spectrum. The shape of the radiation field exciting the 
dust is: 

U,iu)^Ll-^^\u)+Ltiu). (10) 

This dust component is attenuated by the same factor than the diffuse ISRF, P^^^^u), 
in order to reproduce the silicate extinction feature at 9.7 /im. 

5. Our mid-IR spectra exhibit several fine structure ionic lines: [Ar llje.gs^tm, [Ar injs.gg^im, 

[SlV]io.51^m, [NeIl]i2.8lMm, [Nc III] 15.56;,™ , [S III] 18.68/.™, [S IIl]33.48Mm and [Si Il]34.82/.m- We 

fit these lines in order to get a better x^, but we do not use them for our physical 
interpretation. We adopt a Gauss profile whose width is determined by the resolution 
of the spectrograph. The luminosity, Li, of each line is free to vary. 

The total SED is then: 
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The energy conservation implies that: 

{PeLpBR) = -^^burst " (-Pe-^burst) + — {P(,L^) . (12) 

Fig. [5] demonstrates this model, and Fig. [6] illustrates its geometry. 

3.3.2. The Model Results 

We fit the observed UV-to-radio SED, LfW)±/\L°J'W)/2 {i denoting an individual 
wavelength), of each one of the sources in Tabled! with Eq. f|TT]) . minimizing the x^- We 
weight each data point depending on the density of adjacent points: 



, _ ^ min,(A, - A.) ( K^^X^) - L.(A,) \ 
^ Y V AL°b«(A,) J ■ ^ 

Eq. (fT3|) prevents the to be dominated by the parts of the electromagnetic spectrum where 
its sampling by the observations is dense. From this fit, we derive the value of the various free 
parameters. In particular, we are interested in the total PAH-to-gas (and PDR-dust-to-gas) 
mass ratio, ZpAH(dust): 

MpAH(dust) , 1 

^PAH(dust) = ... ^ . . ^ where /i = — . 14 
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Fig. 6. — Illustration of the geometry of the model (Eg. [TTj) . The left panels, from the top to 
the bottom shows the combination of the various SED building blocks from the massive star 
clusters to the observer. The solid lines are the total SED at each step; the dashed lines are 
the SED of the previous step; the part of the SED that has been absorbed is shown in grey. 
The right panel illustrate the path of the photons from the star clusters to the observer. 
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4. THE PAH AND DUST ABUNDANCES IN GALAXIES 

We applied the method described in ^to the sample in Tabled] (Figs. [71 to [Ml) . The 
results are presented in Fig. [25] and Table [3] 

PAH features are not detected in I Zw 18, SBS 0335-052, Mrk 153, NGC 5253 and 
NGC 1399. For these galaxies, we get an upper limit by fitting the short-wavelength part of 
the mid-IR spectrum with the maximum PAH amount allowed by the observational errors. 



I Zw 18: The H I extends far out of the star forming region. Instead of normalizing the 
dust masses by the value of Mhi in Table [T] we consider t hat the H I mass as sociated 
to the star forming region is the one of the object Hl-A ( Ivan Zee et al.l [19981 1. where 
Mh, = 4.4 X 10^ Mf^. 



SBS 0335-052: The 65 /zm flux is the one reported bv lHunt et all (120051 ). The fit of this flux 
gives a far-IR dust temperature slightly colder than what we would obtain by fitting 
only the slope of the Spitzer/YR^ spectrum. However, this gives a conservative solution, 
since we derive only an upper limit on the mass of PAHs, in this galaxy. The H I halo 
extends also far out of the star forming region. We assume that the star forming region 
has the same size than the optical galaxy, i.e. a radius of 1.6 kpc (6"). We derive the 
correspondi ng mass of H I, using t he average column density of iVn i 
reported bv IPustilnik et al.l J200l[ ). We find Mhi = 4.6 x 10^ Mq. 



7.4 X 10^0 cm-2 



VII Z w 40 3: This Spitzer/YR^ spectrum of this galaxy has been studied by [Wu et al. 
(|2006[ ). However, they did not report any PAH detection. The degradation of the 



spectral resolution that we performed on this spectra (see §2.3p increases the signal-to- 
noise ratio significantly. We report a 4cr detection of the 7.7 /im feature, and marginal 
detections of the 6.2 and 8.6 /xm features. This is the lowest metallicty PAH detection 
to date. Similarly to I Zw 18 and SBS 0335-052 the Hi halo extends far out of the 
star forming region. To correct for this effect, we consider that the star forming regio n 
has a size 1' x 0.5', with a column density Ahi = 1.6 x 10^^ cm~^ (IThuan et al.ll2004l ). 
This leads to an effective Mhi = 1.1 x 10^ Mq. This galaxy has only one radio point. 
To remain conservative, we assume that it is free-free dominated. 

Mrk 153: The mid-infrared spectrum of this galaxy show prominent silicate emission. We 
do not have any radio observation of this object, therefore the contribution of Hll 
regions to the total SED is rather uncertain. 



Haro 11: As quoted by [Bergvall et al.[ (|2000[ ). the ratio Mhi/^vb < 0.01 is abnormally low 
in this galaxy. The low value of the H I mass explains the high value of its dust-to-gas 



- 29 - 



mass ratio (Fig. [251) . On the contrary, its PAH-to-dust mass ratio is not peculiar. This 
source, at 92 Mpc, is one of the most distant object in our sample. 

IC 342: This galaxy is located at 10° of Galactic latitude. Hence, it is highly extincted by 
the foreground, which explains the peculiar shape of its optical/near-IR observed SED. 

Tol 89: The mid-IR spectrum of this galaxy samples only its nucleus. That is the reason 
why it underestimates the Spitzer/lRACsfj,m flux. To be conservative, we derive the 
lower limit on the PAH mass by fitting the spectrum, and the upper limit, by fitting 
the broad band. This galaxy has only one radio point. Therefore, we assume that it is 
free-free dominated. 

NGC 1068: This galaxy contains a powerful AGN. In principle our model can not be applied 
to this object, since it does not take into account the contribution of the accretion disc. 
However, we fit this galaxy as if it was a starburst, in order to test the robustness of 
our approach. 

NGC 1399: This object is a cD galaxy. Its interstellar medium is likely very tenuous and 
the mid-IR emission is dominated by the contribution of evolved stars. It has not been 
represented on the upper panel of Fig. [251 since we only know the upper limits on the 
PAH and Hi gas masses. However, it is shown on the lower panel; the higher open 
circle is the ratio between the upper limit on the PAH mass and the lower limit on the 
dust mass. 

Circinus: This galaxy is located at Galactic latitude below 5°. Like IC 342, its optical 
fluxes are very uncertain. 



The top panel of Fig. 1231 shows the variation of the PAH and dust to gas mass ratios 
with the metallicity of the interstellar medium. Each individual galaxy can be seen as a 
snapshot of galaxy evolution, at a given time. First, we note that the Galactic values of 
the dust-to-gas mass ratios are in agreement with the oi ie of the other gala xies, around the 
same metallicity. These Galactic values were derived by lZubko et al.l (120041 ). from the fit of 
the emission and extinction of the diffuse interstellar medium, with further constraints from 
the elemental depletion pattern. Thus, it is a very reliable estimation. This comparison 
confirms that our method does not overlook a significant amount of dust, at least around 
the solar metallicity. Second, the trends of Zp^n and .^dust, with the metallicity are not 
identical. Our sample spreads two orders of magnitude in metallicity. We can see that the 
PAH-to-gas mass ratio rises by five orders of magnitudes, while the dust-to-gas mass ratio, 
by only three. This differential evolution is illustrated in the lower panel of Fig. [251 showing 
the PAH-to-dust mass ratio. This figure is the analog of Fig. [H but instead of considering 
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integrated fluxes, it deals with abundances. The PAH-to-dust mass ratio rises by two orders 
of magnitude, over our sample, while the IRAC8^m/MIPS24^m band ratio varies only by one 
order of magnitude. 
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1 100 10* 10' 

Wavelength [jum] 



Fig. 7. — Fit of the galaxies' SEDs. The black lines are the components of the model (Fig. [5]) 
and the thick grey line is the total observed SED. The circles with error bars are the broad- 
band and spectral observations. The top-right panel shows the detailed fit of the mid-IR 
spectrum. 
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Fig. 8. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 9. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 10.— Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 11.— Fit of the galaxies' SEDs. See Fig. [7] for details. 
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12. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 13.— 



Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 14. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 15.— 



Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 16. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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17. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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18. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 19.— Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 20.— Fit of the galaxies' SEDs. See Fig. [7] for details. 
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21. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 22.— Fit of the galaxies' SEDs. See Fig. [7] for details. 
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23. — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 24 — Fit of the galaxies' SEDs. See Fig. [7] for details. 
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Fig. 25. — Top panel : dust-to-gas mass ratio for PAHs, ^pah? and the dust that gives rise to 
the far-IR emission, Zdust? as a function of metaUicity. Bottom paneh mass ratio of PAH-to- 
dust, as a function of met aUicity. The circle s correspond to galaxies and the open stars to the 
the diffuse Galactic ISM (jZubko et al.ll2004l ). The filled circles are the reliable measurements, 
and the open circles are the ones which are considered uncertain. The grey stripes are the 
±lcr linear correlation, in logarithmic scale. These figures show the two distinct evolutionary 
trends of 2'pah and ^dust with metaUicity. 



Table 3. SED Modeling Results 
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5. A MODEL FOR THE GAS AND DUST EVOLUTION 

We have developed a one-zone single-phase chemical evolution model, to follow the 
abundances and composition of the dust and the metallicty as a function of time, in order 
to interpret the results of our SED modeling. In this section, we will give a brief description 



m 



Dwek 


(1998 


) and 


Dwek et al. 


(2007) 



use this model to interpret the observed evolutionary trend of galaxies' SED with metallicity, 
on a global scale. 



5.1. Metal Enrichment and Gas Evolution 

In the present paper, we consider a closed box model. We consider the delayed injection 
of material by different stellar progenitors, but we assume that the mixing of the elements in 
the ISM is instantaneous. To be consistent with our stellar population synthesis, we adopt 
a Salpeter initial mass function, 0(m) (Eq. [7]), where m is the mass of individual stars. We 
define the average stellar mass: 

rrriu 

{m) = / m(l){m) dm. (15) 

J mi 

The evolution of the gas mass surface density, Egas(t), with the time t, is: 

- ^SFR(t) + / SsFR(t - r(m)) 0(m) dm, (16) 



dt Jmi {m) 

where EsFR(i) is the mass of star formed per unit time and per unit surface area, r(m), the 
lifetime of a star of mass m, and mej(m), its mass of gas returned to the ISM. The first 
term of the right hand side of Eq. (|T6|) is the amount of gas removed by the star formation, 
and the second term, is the delayed injection of the gas (H, He and metals), by the various 
progenitors. The ISM metallicity, is defined as: 

^ISm(^) = ^gas(^) + ^dust(05 (17) 

where Zgas(t) is the metal-to-gas mass ratio, and Zdnst{t), the dust-to-gas mass ratio that 
will be discussed at §5.2[ Its evolution is analog to Eq. ( fT6l) : 

d[Egas(t)^iSM(t)] _ ^ , r% , ,,Yzim 



dt 



Zisuit) SsFR(t) + / SsFR(t - r(m)) 0(m) dm, (1^ 

Jmi {m} 



'mi 

where Yz{m) is the yield of elements heavier than He, by the stars of mass [m, m + dm]. 



Table 3 — Continued 



Name 


Zgas 
[Zq] 


^PAH 


■^dust 


*SF 

[Myr] 


"H /+ 
[cm-3] 


^ burst 

[Lq] 


^PDR 

[Lq] 


M 83 


r, O+0.5 

^••^-0.6 


6.0^4 g X 10-5 


r.llla X lO""* 


5440 


290 0.6 


2.5 X 10^ 


1.7 X 10^" 



Note. — Most of these quantities are defined in ^ and Table |2] Zpp^^l ^^'^ -^dust ^^'^ PAH-to-gas 
and dust-to-gas mass ratios; tgp is the age of the galaxy constrained by the optical and near-IR broad 
bands; nj^ is the average gas density in Hll regions; /+ is the mass fraction of cationic PAHs; ^burst 
recent burst of star formation; L^^'^ is the intrinsic bolometric luminosity omitted by the dust inside the 
neutral phase. The sources arc ordered according to their metallicity. The symbol identifies uncertain 
values. 
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The elemental y i elds, Y z(m), for t he low mass stars (m < 8 Mq) are taken from 
Karakas fc Lattanzid (j2003al Jbl). and from lWoosley fc Weaverl (119951 ) for the high mass stars 
(m > 8 Mq). Furthermore, a prescription for the star form ation rate is req uired, in order 
to solve these equatio ns. This prescrip tion is given by the ISchmidtl (119591 ) law, with the 



coefficients derived by iKennicuttI (119981 ): 



pc 2 



(2.5 ±0.7) X 10 



-10 



^gas(^) 
1 Mq pc- 



1.40±0.15 



(19) 



and is used to replace SsFR(i) in Eq- ( 0^61) . Fig. [26] shows the evolution of the total metallicity. 



and of the reduced gas mass Hgas{t) 
mass surface density. 



where S° „ 

gas 



Sgas(O) is the initial gas 
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Fig. 26. — Evolution of the total metallicity (in units of Zq), and of the reduced gas mass 
(/igas), for two different initial surface gas mass densities, S° . 



5.2. Dust Formation and Destruction 



The evolution of the mass surface density, of a given dust specie, is the sum of 

three contributions: (i) the rate of dust destruction by star formation, (ii) the rate of dust 
condensation in stellar progenitors, (iii) the rate of dust destruction in the ISM, by SN blast 
waves: 



dSi(t) 



-Z,{t) EsFR(t) + 



mi 



Yi{m) 
(m) 



^SFR 



{t — T{m)) (j){'m) dm 



Tdust(^) 



(20) 



where Yi{m) is the yield of the considered dust specie by stars of mass [m,m + dm], and 
'7"dust(^)5 the dust lifetime in the ISM. 

T he dus t yields are derived from the elemental stellar yields described at §5.11 following 

Dwekl (119981 ). For low mass stars (m < 8 Mq), the dust yields depend on the value of the 
C/0 ratio. We assume that the ejecta is microscopically mixed, so that all the excess carbon 
is locked-up in dust, if C>0. If C<0, then we combine all the available Fe, Si, Mg, Ca and 
Ti, with one O atom to produce silicate dust, and titanium oxydes. For high mass stars 
(m > 8 Mq), we assume that the ejecta is only macroscopically mixed, so that both carbon 
and oxygen rich dust can condense. We assume a condensation efficiency of unity for all 
dust species. 



The life time of dust, rdust(^) is directly related to the SN ll rate (IDwek fc Scald Il980l : 
McKeelll986[ ): 



Sgas(^) 



(21) 



RsN ii(^) ("Iism) ' 

where -RsNn(^) is the SN ll rate per unit area, and (misM) is the average effective mass of 
gas, swept up by a single SN remnant, where the dust is returned b ack to the gas phase, by 
either thermal sputtering, or grain-grain collision (j Jones et al.lll996l ). The value of (ttiism) is 
largely unknown. We explore its effects, by varying it from (misM) = Mq (no destruction), 
to (misM) = 300 Mq (typical destruction). The latter value is typical of our Galaxy. Indeed, 
if M5,^^ ^ 5 X 10^ Mq is the total mass of gas in our Galaxy, R^^^ ~ 1/3 yr~^, its a verage 
5 X 10^ yr, the mean lifetime of an ISM dust particle (jJonesll2004l ). then: 



gas 



SN rate, and r^f^ ~ 



("^ism) 



gas 

f?Gal _Gal 
-"-SN Must 



300 Mq. 



(22) 



An additional destruction mechanism that apllies only to PAHs and very small grains 
is their photoevaporation in intense radiation fields. This destruction mechanism primarily 
affects the abundance of PAHs, and is taken into account in the SED model that calculates 
the dust abundances, by assuming that PAHs do not survive in H ll regions. 
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Fig. [27] shows the evolution of the total dust content formed by massive stars, and the 
carbon dust formed by AGB stars, for various initial conditions, and destruction efficiencies. 
The SN II dust evolves almost linearly with the metallicity, in absence of destruction, because 
the metal enrichment is dominated by massive stars. In contrast, the AGB carbon dust starts 
rising when the metallicity of the ISM is around 1/20 Zq. This value corresponds to a time 
of ~ 100 Myr (Fig. [26l) . which the lifetime of the most massive AGB stars. Hence, the carbon 
dust produced by AGB stars is injected into the ISM, with a delay whi ch correspon ds to 
the lifetime of the st a rs. T his evolutionary trend was previously noted by iDwekl (119981 ) and 



Morgan fc Edmundsl (120031 ). The change in the slope of the AGB carbon dust, around 1 Zq, 
is simply due to the fact that AGB stars of lifetime longer than ~ 1 Gyr are oxygen rich. 
The dust destruction effects the evolution for Z > 0.1 Zq. 




Fig. 27. — Evolution of the dust content with the metallicity of the ISM. We adopt (misivi) = 
300 Mq, for the curves with grain destruction. 
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5.3. Examining the Consistency Between the Stellar Populations and the 

Elemental Evolution 

For consistency, we used the same prescription for the evolution of the star formation 
rate in both the population synthesis (§!]) and the elemental evolution modeling ( §5.ip . 
However, we added a late burst of star formation to calculate the galactic SED, which will 
only have a limited effect on the final elemental abundances of the galaxy. The two models 
differ in the stellar yields, with the latter using the more recent yields for AGB stars. Fig. [28] 
examines the effect of this different yields by plotting the galactic age versus metallicity 
relation for our sample of galaxies. The model calculations, shown as a grey stripe, are 
in general good agreement with the data. However, for sub-solar metallicities (< 0.3 Zq), 
several galaxies appear to have an older age, from the stellar point of view, than what would 
be inferred from their metal enrichment. This discrepancy is probably the consequence of 
assuming a smooth star formation history. The differences will only manifest themselves as 
a moderate change in the stellar spectra, and will have no affect on the conclusions of the 
paper. 



5.4. PAHs and the Delayed Injection of AGB Carbon Dust 



From an observational point of view, PAHs are believed to form in the circumstellar 
envelopes of carbon rich AGB stars, and to be subsequently ejected into the ISM through 
stellar winds. The paucity of UV photons prevents the direct detection of these P AHs, during 



their formation process in the post-AGB phase of the evolution of these stars flHony et al. 



2001: Boersma et a. 



20061). However^ PAHs are abundantly observed during the later plan- 



etary nebula phase 
formed in the envelopes of AGB stars. 



Hony et al.l (120011 ). In what follows, we will assume that PAHs are only 



Fig. [22] shows the comparison between the dust and PAH abundances derived from 
the observations (Fig. [25]) . and the ones produced by the evolution model (Fig. [27j) . The 
agreement between the PAH-to-gas mass ratio, and the carbon dust production by AGB 
stars is very good, with the galaxies He 2-10 and IC 342 as the only outliers. He 2-10 is a 
solar metallicity blue compact dwarf galaxy. Its SED is similar to lower met allicity galaxies 
([Galliano et al.ll2005h . ^lowever its metallicity is very uncertain; for example IVacca &: Contil 
( I1992I ) quoted 1/6 Zq. This uncertainty may reside in the fact that this galaxy has two 



cores, result ing from a merger. The properties of these nuclei are differen t in terms of dust 
absorption ( Phillips et al. 1984), m olecular gas content ( Baas et al. Il994), compact sourc e 
distribution (jCabanac et al.ll2005l ). and mid-lR spectrum (IMartin-Hernandez et al.ll2006l ). 
Thus, our global approach may not apply to this object. On the contrary the PAH-to-dust 
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Fig. 28. — Consistency between the stellar populations and the elemental evolution. The 
error bars are the age of the galaxies derived from population synthesis modeling of the optical 
and near-IR broad-band observations (Table [3]), as a function of the observed metallicity of 
the gas (Tabled]). The grey stripe shows the range of values from the elemental evolution 
model (Fig. [26]). 
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Fig. 29. — Comparison between the metallicity trends of the PAH abundance derived from 
the observed SED and those derived from the chemical evolution model. The figure highlights 
the different evolutionary trend of SN ll- and AGB-condensed dust. 
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mass ratio of IC 342 is consistent with the one of other galaxies with the same metaUicity 
(Fig. [25]) . Only its values of 2'pah and ^dust are systematically shifted. As mentioned in 
§1] this galaxy suffers a lot of extinction and confusion since it is located at low Galactic 
latitude. This could alter the estimate of the absolute dust-to-gas mass ratios of this object. 
We conclude that the striking evolution of the PAH content in galaxies with metallicity can 
naturally be explained by the delayed injection of carbon dust into the ISM, by AGB stars. 

In principle, PAHs can form by other processes, for example, by the hydrogenation of 
small carbon grains in the ISM. Initially, their abundance will then follow the evolution- 
ary trend of SN-condensed carbon dust. However, their absence in low metallicity systems 
suggests that they are efficie n tly de str oyed, presumably by sho cks or UV radiation, as en- 
visioned in by iMadden et al.l (120061 ) or lO'Halloran et al.l (120061 ). At later times, the rate of 
PAH production by these interstellar processes will follow the evolutionary trend of AGB 
stars, since they will be the major source of carbon dust in the ISM. 

The global trend of PAH abundance with metallicity does not preclude the possibility of 
local variations of PAH abundance in indivi dual galaxies. For e xample, the SMC which has 
on the average a metallicity that is 1/6 Zrr, (iDufour et al.lll982l ) has at least one region, the 
molecular cloud SMC Bl#l (Reach et al. 200ol ). with a PAH-to-dust ratio that is comparable 



to the Galactic value (ILi fc Draindl20021). Even in our Galaxv. there are regions with "super 



Galactic" PAH abundances (e.g. lRidderstad et al.ll2006l ). Such local abundance variations are 



a natural consequences of processes, such as mixing of stellar ejecta and cycling between the 
ISM phases, that when globally averaged, will follow the general trend of PAH abundances 
with metallicity depicted in Fig. | 



5.5. The Evolutionary Trend of SN-Condensed Dust 

The dust content is in good agreement with the dust production by SN ll, down to 
0.1 Zq fFi g. 12^. Below this value, th e model systematically overestimates the observed 



dust content. iLisenfeld fc Ferraral (119981 ) found a similar deficiency, from IRAS observations 
of blue compact dwarf and dwarf irregular galaxies. There may be several reasons for this 
discrepancy. 

First, we may have overestimated the gas mass, in the lowest metallicity sources. Indeed, 
the distribution of Hi of these galaxies extends farther out of the star forming region, as 
discussed in §H We have attempted to correct this effect, for I Zw 18, SBS 0335-052 and 
VII Zw 403, by considering only the gas content associated with the optical galaxy. However, 
we could not correct for the gas located out of the star forming region and along the line of 
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sight. 



Second, we may have underest i mate d the dust mass by overlooking a cold dust compo- 



nent. Indeed, iGalliano et al.l (120031 . |2005| ) showed that the millimetre excess observed in the 
SEDs of NGC 1140, NGC 1569, II Zw 40, and He 2-10, could be consistently explained by 
the presence of very cold dust, accounting for 40 to 80% of the total dust mass. On Fig. [291 
the dust-to-gas mass ratio of these four galaxies is indeed, below the SN ll production rate 
by a factor of ~ 2. Assuming that the high dumpiness of the ISM is a general property of 
low-metallicity systems, and that the filling factor and/or contrast density of these clumps 
rises when the metallicity drops, we have a natural explanation for this deviation. We can 
not address this issue, because of the lack of submillimetre data for these very low-metallicity 
galaxies. 

Third, the discrepancy between the predicted and observed SN ll dust at low metallicity 
could be due to several parameters or assumptions used in the chemical evolution modeling. 
(1) The IMF could play an important role in the absolute value of the dust production by 
SN II. It would indicate that the slope of the IMF is metallicity-dependent, and that the 
contribution of massive stars is lower, at very low Zisu, which seems to be unlikely, both 
from a theoretical point of view (star formation), and an observational point of view (star 
counts). (2) A much more likely explanation could come from the fact that we assumed 
that the condensation into dust of the elements ejected by the SN ll and their mixing in the 
ISM is instantaneous. If the major part of the dust was to condense into the ISM, and not 
directly into the SN ll-ejecta, then the dust formation would be delayed after the death of 
the massive stars. (3) Finally, the dust production rates, computed from our dust evolution 
model, implicitely assume that the star formation rate of the galaxy is smooth and that the 
mixing is instantaneous. This hypothesis is ce rtainly valid fo r evolv ed systems, but could be 



wrong for very young objects. For example, iLegrand et al.l (120001 ) suggested that the star 
formation history of I Zw 18 is not continuous. 



6. CONCLUSION AND SUMMARY 

The weakness of the mid-IR aromatic features in low-metallicity environments has been 
traditionally interpreted as the consequence of the increased selective destruction efficiency 
of the PAHs in these environments. In this paper, we presented a new interpretation for 
the observed correlation of the intensity of the mid-IR emission from PAHs with galaxies' 
metallicity. In our model, this trend is a manifestation of the evolution of the abundance 
of interstellar carbon that formed in AGB stars with time. A fraction of the carbon dust 
formed in AGB stars in the form of PAH macromolecules offering a natural explanation of 
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the evolutionary trend of PAH abundance with galactic metallicity (or time). 

To ascertain this trend, we first converted the trend of observed mid-IR fluxes to PAH 
and dust abundances, by modeling the UV-to-radio SEDs in a sample of 35 nearby galaxies, 
with metallicities ranging from 1/50 to ~ 3 Zq. Our models represent the most detailed 
decomposition of the dust emission into its gas phase components using a wide range of 
astrophyical constraints, including: the free-free and mid-IR emissions to constrain the gas 
and dust radiation from H II regions; the far-IR and optical emission to constrain the ISRF 
that heats the dust in PDRs. From these models, we determined the abundances of the 
PAHs and other dust species associated to the neutral phase of the ISM. We then used a 
chemical evolution model to calculate the abundances of SN- and AGB-condensed dust as 
a function of time or metallicity. The model takes into account the delayed recycling of the 
ejecta from low-mass stars caused by their finite main-sequence lifetime. 

The main conclusions of this paper are the foUowings. 

1. From the decomposition of the IR emission into its various emission components and 
dust species, we discovered two distinct evolutionary trends for PAHs and other dust 
species in galaxies. 

2. Comparison of these trends with the chemical evolution models showed that the abun- 
dance of PAHs follows that of the carbon dust from AGB stars. The remaining dust 
species follows the evolutionary trend of SN-condensed dust. 

3. The delayed injection of PAHs into the ISM provides a natural explanation for the 
paucity of these large molecules in low metallicity systems. The subsequent rise in the 
PAH-to-gas mass ratio with metallicity is then a natural consequence of the increasing 
contribution of AGB stars to the chemical enrichment of the ISM as they evolve off 
the main sequence. The trend of the other dust species is a natural consequence of the 
evolution of SN-condensed dust which is instantaneously injected into the ISM after 
the birth of the progenitor star. 

4. While the model is very successful in observing the gross general trend, there are some 
systematic deviations showing that dust-to-gas mass ratios inferred from observations 
fall below the calculated value for SN-condensed dust in the lowest metallicity galaxies. 
This discrepancy may be due to one or more of the following: (i) a.n overestimate of 
the Hi gas mass used to derive these dust-to-gas mass ratios; (ii) an underestimate of 
the dust mass due to the possible presence of a cold dust component; and (in) a more 
complex star formation history than used in the model calculations. 
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5. As an aside, in our analysis of the mid-IR spectra of nearby galaxies, we detected 
the 7.7 /im aromatic feature at the 4a level in Vll Zw 403, a dwarf galaxy with 
Z ^ 1/20 Zq. This is the lowest metallicity galaxy for which PAH emission has been 
detected to date. 

The success of our chemical evolution model, in reproducing the trend of PAH abun- 
dances with metallicity, strongly suggests the importance of stellar evolutionary effects in 
determining the abundances and composition of dust in galaxies. These will have impor- 
tant consequences for determining the opacity of galaxies and their reradiated thermal IR 
emission. Chemical evolution models for dust must therefore be an integral part of popula- 
tion synthesis models, providing a self-consistent link between the stellar and dust emission 
components of the SED of galaxies. 
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